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Abstract 

We calculate the lifetime of the collective isospin rotating solutions which 
have been found recently in the case a quantum N-component meson field 
with exact 0(N) symmetry. For this purpose we take into account the small 



\j^ \ breaking of the 0(N) symmetry associated to the non vanishing mass of the 

C^ I pion. This term induces a coupling between collective rotations and intrinsic 

(^ ■ meson excitations. We evaluate the associated damping time in the framework 



of linear response theory. We find damping times of the order of 100 fm/c, 
i.e. substantially longer than reaction times. 
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I. INTRODUCTION 

The forthcoming generation of ultrarelativistic nuclear colhders at Brookhaven and 
CERN opens the perspective of exploring the dynamics of matter at very high energy density. 
Among the most central issues in this field is the production and study of the quark-gluon 
plasma. Another attractive question is the physics of high density meson fields out equi- 
librium which could evolve on a time scale of the order of a fm/c or less. Such systems 
are of interest for the understanding of the field theoretical problems involved in the first 
stages of the early universe. They may also exhibit collective phenomena closely connected 
to fundamental symmetry properties such as chiral symmetry. The formation, the stability 
and the damping time of such collective modes are important, yet unanswered questions, 
although some progress in this direction has been carried out recently [|T]-§] • 

Collective modes in a high density pion field have been studied in the context of clas- 
sical evolution equations by Anselm and Ryskin [Q and also by Blaizot and Krzywicki |^ 
who focused on the issue of Lorentz invariant boundary conditions. These time-dependent 
solutions of classical meson fields attracted wider attention from the point of view of non- 
equilibrium dynamics of the chiral phase transition §,0]. Recently quantum fluctuations in 
the description of these modes were introduced 0,3. Particular analytic solutions of the 
mean field evolution equations corresponding to quantum rotations in isospin space were 
constructed, both in the zero temperature case as well as in the finite temperature case 

In the limit of exact chiral symmetry these rotations decouple from intrinsic meson exci- 
tations. However when an explicit symmetry breaking term, which generates non zero pion 
mass, is introduced, then this term will act as an anisotropic (in isospin space) external po- 
tential for the isospin rotation which causes coupling of the collective rotation and intrinsic 
mesonic excitations of the system and this coupling in turn results in a dissipation of the 
collective excitation energy to mesonic excitation energy. This effect is very similar to the 
Landau damping of collective oscillations in collisionless plasma. Indeed, as we shall see 
below, the symmetry breaking term is viewed in the isospin rotating frame as a periodic 
time-dependent external perturbation, the linear response of the system to such an external 
oscillatory source contains collisionless dissipation, a characteristic feature of the Landau 
damping mechanism generalized to other systems. As a result the collective "angular mo- 
mentum" for the isospin rotation diminishes. Since the collective isospin angular momentum 
may be proportional to the "moment of inertia" of the isospin rotation (the square of the 
condensate) times the corresponding angular frequency, both the angular frequency and the 
radius of the chiral circle, which are related each other in the gap equation, may become 
small. 

It is the purpose of the present work to evaluate the corresponding damping time. In 
order to perform our calculations we will assume that the coupling between collective ro- 
tations and meson excitations via the symmetry breaking term is small. This is the case 
when the kinetic energy density of collective rotations is large compared to the change in 
the potential energy density induced by symmetry breaking term. This can be seen to be 
the case when the frequency u of rotational motion is large compared to twice the pion mass 
(see below). 

In the following we will investigate properties of the damping time such as its depen- 



dence on excitation energy density, angular frequency and temperature. Another question 
of interest is to work out the spectrum of emitted particles. Experimental issues related 
to the previous questions are the observability of classical meson field as suggested first by 
Anselm and Ryskin [^]. In particular it is obviously important to have a reliable estimate 
of the lifetime of the collective modes of the pion field in order to be able to check that it is 
longer than the reaction time. Also it is of great interest to have estimates of the lifetimes of 
various modes in order to find out whether damping could act as a filter to select the most 
collective states as the nuclear collision proceeds. This information is crucial to discuss the 
formation of collective modes. 

In section 2 we review the isospin rotating solutions of the quantum mean field equations 
which were obtained in reference 0, hereafter referred to as I. These solutions correspond to 
static solutions in the isospin rotating frame in the case of exact chiral symmetry. In section 
3 we consider the case where there is an explicit symmetry breaking term. We show that this 
term generates a periodic external source in the rotating frame and we construct an explicit 
expression of the source. In section 4 we describe the general formalism of linear response 
theory in the rotating frame which is suited to deal with the case of a small symmetry 
breaking. As a warm up exercise we deal in section 5 and 6 with the simple case of the 
bare response in which residual interactions are neglected. The realistic situation where 
interactions are included is discussed in section 7. The expression of the damping time is 
derived in section 8 and a discussion of its physical content is presented in section 9. In 
Appendix A, we investigate the pionic excitation in the chiral limit in the linearized mean 
field approach and show that our approach is consistent with the Goldstone theorem. 

II. REVIEW OF ROTATING SOLUTIONS OF THE MEAN FIELD EQUATIONS 

Time-dependent solutions of the classical equations of motion of the non-linear sigma 
model have been constructed and discussed by Anselm and Ryskin [Q and others p . These 
solutions describe field rotations in its internal symmetry space which may be considered 
to appear as a result of the spontaneously broken internal symmetry, similar to nuclear 
deformation and rotation. Here we review a quantum extension of the Anselm-Ryskin 
solutions in the framework of the mean-field description as obtained in [^ and in I. 

For N-component scalar fields (pa{x) with a bare mass tUq and a coupling constant A, the 
evolution of the mean field if is governed by [Q 

{□ + TTig + -(p"^ + -traceS'(x, x) + -Six, x)}(f = 0, (1) 

DO 6 

where the trace runs over the index a and the NxN matrix operator S is the Feynman 
propagator, iS~^ = D + m^(x) + ie, with the mass matrix given by 

m^(a;) = 1711-1 — <^^(x) H — trace5'(x, x) H — (^(x) x i^{x) H — S{x, x), (2) 

6 6 3 3 

where the symbol <^(x) x ip{x) denotes the NxN matrix whose matrix (a, 6) element is 
i^a{x) X i^b{x). To construct solutions to the above coupled equations, we adopt the following 
Ansatz: 



(p{x) = U{x)ipr = exp{i{ujt — q ■ x)Ty}ipr, (3) 

and 

S{x,y) = U{x)Srix,y)U\y) , 

where (pr is the amphtude of the mean field in the "rotating frame", Ty, the second Pauh 
matrix, is the generator of the rotation in the subspace of 0(N), 

■•\ 

:: . (4) 

■■) 

and uj and q are arbitrary parameters to characterize space-time dependence of the rotation. 
The propagator Sr in the rotating frame takes the form of 

Sr{x,y) = ^g^ ^ ^q^ry){^^ + ^q,ry) + M^ + is' ^^^ 





/O 


—i 





^y = 


i 












where 



M^ = U{x)m\x)U^{x). (6) 



In the rotating frame, the condensate ipr and the mass matrix M^ satisfy the mean field 
equation (|I]) and the gap equation (0) with S{x,x) and (f{x) replaced by Sr{x,x) and i^r, 
respectively, and one can show that M^ becomes independent of space and time. These two 
equations form a closed set to determine the values of M and (pr self-consistently. 

In the limit q=0 one recovers the self-consistent equations specifying the properties of the 
vacuum state, namely the associated values of the mass matrix m and of the condensate y^o- 
These equations are reproduced in Appendix A (equations ( |A12| ) - ( [A15| )). The mass matrix 



m is diagonal in the vacuum state, with A^ — 1 meson masses being equal to a value jj, (which 
corresponds in the sigma model A^ = 4 to the pion mass), the last one being equal to M 
(the sigma meson mass in the sigma model). It is worthwhile noting that even in the chiral 
limit (no symmetry breaking term) the pion mass in the mean field approximation does 
not vanish (see Appendix A). However the existence of rotating solutions guarantees that 
the linearized mean-field evolution equations do have solutions with zero frequencies. As a 
result the pion masses in this more elaborate framework do vanish as shown by an explicit 
calculation described in Appendix A. This implies that the Goldstone theorem is recovered 
only at the level of the linearized mean-field evolution equations. This is in agreement with 
the discussions already given by other authors [p!0H l^. 

The solutions we have just reviewed concern only the zero temperature case, i.e. pure 
states. Generalization of these solutions to finite temperature was reported in I. In this case 
one needs to consider statistical mixtures. These are most conveniently handled by introduc- 
ing the corresponding density operator formalism |14,1^. In this formalism the mean field 



approximation is formulated simply by stating that the density operator D describing the 
system is given at each time t by the exponential of a quadratic form in the field operators. 
In this case the expectation value of an arbitrary operator 



(O) = TracejDO} (7) 

can be expressed in terms of the generalized density matrix Ai defined by 

Mix,y;t)+^-=( -'[^^^'^l^f^f f'^[''\^.'^y\\ ) . (8) 

^ '^' ^ 2 V {na{x)7rt{y)) i{'Ka{x)^b{y)) ) ^ ' 

In this equation 77^(0;) is the momentum conjugate to '^h{p^\ We have used the notation 

(paix) = ipa{x) - {(Pa{x)). (9) 

Expectation values of field operators are defined by equation (^. The mean field evolution 
equation for the generalized density matrix is obtained in this formalism by stating that the 
trace of the Liouville equation for D multiplied by an arbitrary quadratic form in the field 
operators vanishes, i.e. 

Trace {{iD - [H, D])0aix)0biy)} = 0, (10) 

where H is the hamiltonian of the 0(N) model whose density reads 

^0(N)(a3) = -na{x)'Ka{x) + -V Lpa{x)V LPa{x) + -mlipa{x)Lpa{x) + —{Lpa{x)(pa{x)f . 

By calculating the commutators in this formula one obtains the following evolution equation 
for the generalized density matrix 



iM{t) = [n , Mit)] , (11) 



where 



^-(? J)' (12) 

r"'= l-A + ml + -^^ + -{0,0,)]6ab + -{Va^b+{0a0b)) , (13) 

and 1 is an A^ X A^ unit matrix in the isospin space. The detailed derivation of this equation 
is given in I. 

Static solutions of these equations corresponding to thermal equilibrium at a given tem- 
perature T are obtained by writing that the generalized density matrix Ai commutes with 
the mean field hamiltonian 7i. In this case one finds a Boltzmann-Gibbs type distribution 
for the density operator 

D = ie-^^ , Z = Trace e"^^ , (14) 

Zj 

where /3 = Xjk^T . In the previous equation W is the Hartree-Bogohubov Hamiltonian 

/r 1 1 Tfi 

d?'x\-T^a{x)T^a{x) + -^ (f) a{x)S] (f) ^{x) + -^(pa{x)(pa{x) 

■^-r^{^a{x)^a{x))(pi{x)(p},{x) + -l^i^ a{x)^i,{x)) (p a{x)(pi>{x) \ . (15) 



In order to build rotating solutions at finite temperature it is convenient to write the 
equation of motion for the generalized density matrix in an isospin rotating frame charac- 
terized by a frequency uj and a wave vector q as defined by equation (^. This is done again 
by gauge transforming the density matrix: 

M{x, y; t) = U{x, t)Mr{x, y; t)U\y, t) . (16) 

Then the Liouville equation is transformed to 

iMr = [Hr , Mr], (17) 

where the mean field hamiltonian in the isospin rotating frame is given by 

n. = uH.)(-i^ + H)u(.) = {^;^ j.J. (18) 

rf = I -(V - iqTy)'^ +ml + -ifl + -{^c^c)r j ^ah + ^(<^r,a<^r,6 + {0a0b)r) ■ (19) 

In the above expression, the subscript r indicates that the variable is taken in the isospin 
rotating frame. In the following, we omit these subscripts unless it is very confusing. 

We construct the isospin rotating solutions at finite temperature by imposing that the 
density matrix takes a static, thermally equilibrated form in the rotating frame. In this the 
case Tir and A4r commute and therefore they can be diagonalized simultaneously: 

Here, the static and thermally equilibrated case is described with the subscript 0. As is 
expected, {u*, —v*) is also an eigenvector with an negative eigenvalue —E and — /. Here let 
us introduce a diagonal basis set {\a, k, r)} for TYq and A4q in momentum space, where a and 
k represent isospin index and momentum, respectively, and r = + or 1 (— or 2) represents 
the state with a positive (negative) eigenvalue. Then we obtain 

no\a,k,T) = E:{k)\a,k,r) ; E^{k) = ±E,ik) , 
Mo\a,k,T)=ra{k)\a,k,T); /±(fc) = ± (n„(fc) + ^) , (21) 

where Ea{k) can be directly calculated from Tio and is given in I and na{k) = 
l/{exp{PEa{k)) — 1) is a bose distribution function. 

The density operator which gives a thermal distribution in an isospin rotating frame can 
be shown to be of the form 

D{uj,q) = Le-f^^i-^'^) , Z = Trace e"^^ , (22) 

Zj 

which is the same as before except for the presence of a modified Hartree-Bogohubov Hamil- 
tonian 

2 

W{yj, q) = W-ujIab + q-Iab + Y J d^x0c{x)ipc{x) . (23) 



In this equation the operator I is the generator of isospin rotations 

d^x[ipa{x)Vipb{x) - (pb{x)V(pa{x)] . 



lab — 
lab 



(24) 



We choose a = 1 and 6 = 2 for the rotation in the isospin space in accord with the ansatz @. 
In what follows we will denote by A4Q{uj,q,T) the associated generalized density matrix. 

III. CLASSICAL ISOSPIN ROTATIONS WITH SMALL SYMMETRY BREAKING 



In this section we consider the case where a small explicit symmetry breaking term is 
added to the Hamiltonian density, i.e. 



n{x) = no{N){x) - Cipa{x)Sal ■ 



(25) 



We wish to work out the effect of the second term on our rotating solutions described by 
Aio{uj,q,T) in the rotating frame. In a first step let us investigate the change it produces 
on the isospin rotating solutions obtained in the context of classical dynamics by Anselm 
and Ryskin M . The classical equation of motion is derived as 



d^d'' + ml + ^ipA q> = c\ 



(26) 



We are interested in the isospin rotating solution. In rotating frame, the expectation value 
(condensate) of the field is written as 



<^ = e*9-"«((^(°) + 5¥p(x)) 



(27) 



with q = {lj, q). Here we introduce the condensate ip^, which does not depend on space and 
time, pointing in the direction of the first isospin : 



V 



(0) 







Vo; 



(28) 



The condensate y^o satisfies a gap equation without c in the chiral symmetry broken phase : 

A 



- g' + ^0 + Q <^o = 



(29) 



The quantity 5ip{x) appears due to non- vanishing c, that is, the existence of the explicit 
chiral symmetry breaking term. Then the equation of motion for the condensate is reduced 
to the following equation of motion for 5ip{x) : 



[(5^ + iq>.Ty){d^' + iq''Ty)]ij + {ml + -ipl)5i, + -#<^f \5^i{^) = c 



( cos qx \ 

singx 



V ) 



(30) 



This equation can be solved up to the order of c : 



/( 



5ip{x) 



c 

'a 



[ml 



4g^) cosgx^ 
4g^) sin qx 




V 







/ 



A 2 2 o 2/ 2 I 2 



(31) 



(32) 



where m^ = m^ + (A/6)v9q and m^ = mg + (A/2)(pg. 5ip{x) depends on both time and space 
in the isospin rotating frame. This solution produces a term which plays a role of an external 
source for collective isospin rotation of quantum meson fields. 

IV. LINEARIZED MEAN-FIELD DYNAMICS IN THE ISOSPIN ROTATING 

FRAME 



In this section we examine the effect of the symmetry breaking term on collective isospin 
rotations in the mean field approximation. For this purpose we consider mean field dynamics 
in the presence of the source term derived in the previous section. The boundary condition 
we consider is the following. We assume that at time t=0 a collective rotating state has been 
created, i.e., we take the density matrix in the rotating frame at time t=0 to be M.q{uj^ q, T). 
We also assume that the symmetry breaking term is small enough to allow the use of linear 
response theory. Writing the generalized density matrix Ai{t) as 



M{t) = Mo + 6M{t) , 
and retaining only first order terms in 6Ai in the evolution equation, one finds 



i6M = [Ho , SM] + [6n , Mo] + [SHe^t , M^ 



where 



and 



sn 



/ 

Ur 



SK 



6r 



ab 



ext 



A 







- {6{ipc0c)) Sab + l^H0a0b) , 
6 



6T:^, = ^VoSVlSab + ^(V^i°)5V^6 + ^V^aV^^) 



'{-) 



'{+)/ 



Kb'iQ)e'''-''e^'^' + X:,t>{q)e 



-iq-x iidt 



while 



X 



il\^) 



QA 



.m^ 



V)(2 + r,),,±2(m2-V)(r, 



x/ab 



(33) 



(34) 



(35) 



(36) 



(37) 



Here r^ and r^ are N x N matrices in the isospin space defined as 



Tc: 



fO 1 •■■^ 

1 •■ 

■■ 

Vi ■ ■ •••/ 



rz 



/I ■■■^^ 

-1 ■■ 

■■ 

Vi ■ ■ •■•/ 



(38) 



It should be noted that dHext appears due to the shift of the condensate (mean field) as- 
sociated with non-vanishing c. We regard this term depending on both time and space as 
a classical source term associated with the existence of the explicit chiral symmetry break- 
ing term in the isospin rotating frame. We will investigate the response for this source in 
the succeeding sections. Also, 6TC appears because of the coupling of the intrinsic meson 
excitations and the classical rotating mean field in the isospin space. 

In the basis which diagonalizes both Tio and A^o the linearized equation for 6Ai is simply 
expressed as 



i{a, k, T\SM\b, k', t') = {El{k) - El\k')){a, k, T\SM\b, k', t') 

Hfb'ik') - /:(fc))(a, k, T\Sn + Sn,^,\b, k', r') . 



(39) 



Here the relation between original matrix elements in (|35D and the new matrix elements in 
the diagonal basis is given as 



{a,k,T\6n(^^t)\b,k',T') 



where 



U{k) 



ar.lj 



^-'^'yisTjk,k') o)-^(^') 



<(fc) ui*{k) 
viik) -vi*{k) 



arfir' 



(40) 



(41) 



V. CONSTRUCTION OF THE BARE RESPONSE 

A usual step in the calculation of the full response function to an external field is the 
construction of the so-called bare response function. This function corresponds to the case 
of a simplified evolution equation in which the mean field hamiltonian 7i(t) is approximated 
at each time t by the equilibrium mean- field hamiltonian TYq, which gives 



iM = [Ho ,M]. 



(42) 



For this evolution equation it is straightforward to calculate the response of the system 
(assumed to be in the state A4o at time t = — oo) to a small external field of the form 



6T 



ext 



Xabe 



ici-:x. —iujt 



(43) 



where the frequency uj is assumed to contain a vanishingly small imaginary part to guarantee 
an adiabatic switching of the field. Because of translational invariance in time and space, 
the change in the density matrix at time t is of the form 



6{0a{x)^h{x)) = Yabe'^^^'e 



-iujt 



(44) 



When the evolution preserves isospin, then the matrix Y is also proportional to X. The bare 
response function is defined as the limit of the proportionality coefficient when the external 
field strength tends to zero 

ni?(^,q)=limy/X (45) 

A — >\j 

In this formula we have made explicit the dependence of the bare response function 11*^°^ on 
the frequency cj, on the momentum q and on the isospin indices a and h. 

To obtain an explicit form of the bare response we take matrix elements of the evolution 
equation in the basis which diagonalizes both T^q and A^o- We obtain for the changes in the 
elements of the density matrix 

{a,k,r\5M\h,k\r') = -Jl^f}-J^^}^^{a,k,r\5H.^,%k\r'), (46) 

where the matrix elements of the external hamiltonian read 

'{k)\ 



{a,k,T\5no^t\b,k,T') = X, 



Im 



uT' 



l'^^[^-^){<{k%ut{k'))\j{k'-k + q)e-'-' . (47) 



The evaluation of these elements requires the knowledge of the mode functions ^^(fc), defined 
by the eigenvalue equation (^) and obtained in analytical form in I. Also, we can derive 
f^(fc) by usingof M^(A;) and (^) : v^^{k) = Ea{k)u''^{k)—uJ2b{'Ty)abUbi^)- For more simplicity 
we use in what follows the approximate form which corresponds to the limit g = : 



nUk)^-6,i E,{k) = ^e + M^,. (48) 

/2K{k) 



We also use the approximate values of the masses corresponding to g = 0, namely, Mi = rria 
and M2 = M3 = ■ ■ ■ = Mn = m^, when we estimate the damping time of the collective 
isospin rotation. 

By combining these equations we find that the change in the density distribution reads 



5{ipaix)(fbix)) = Xim e'"^"^e 



iq-^-ic.t I d^k na{k) + nb{k - q) + I 



(27r)3 uj-{Ea{k) + E,{k-q)) 

xui{k)ui^*\k)ut^*\k - q)uf{k - q) , (49) 



where the repeated indices /, /', m and m' are summed over. By using the approximate 
expression of the mode functions we obtain the following expression for the bare response 
function 

r(o)/, . „\ _ f d^k 1 na{k) + nb{k - q) + I 



Cau;,q) 



(27r)3 AEa{k)E,{k -q) uj - {E^ik) + E^{k - q)) ' 

VI. IMPLEMENTATION OF BOUNDARY CONDITIONS 

In this section we consider the response of the generalized density matrix for the exter- 



nal source term (57iext defined by equations (^, ^3) with the boundary conditions already 

10 



described in section 4. The interaction term STi induced by the couphng between the 
collective-rotating condensate and intrinsic meson excitations is still omitted here and will 
be taken into account in the next section. We split 5fA into 5AA^~^ +5M.^^'^ in which (JA^*-^-' 
corresponds to the response for STi^^l oc e^*'^*. Since 6Ai.^^^ oc e^*'^*, we can easily obtain 
the solution for (ROl) without 57i as 



,fc,r|<5^W|6,fc',r') = ^\^^^ ^^^l^l {a,k,T\5ni^^\b,k',T'), (50) 



-[rffc))("-(^')'"-(^'))] /(fc' - fc T q)e^^-* . (51) 



where, from 



Here the factor 6{k' — k^q) reflects momentum conservation arising from the matrix element 
{k\e^''i-''\k'). 

However this particular solution of the linearized evolution equations does not correspond 
to the desired boundary conditions 6M. = at time t = 0. As discussed already in section 
4 this is adequate to describe a situation in which a collective rotating state is created at 
time t =0. In order to satisfy this boundary condition we add a particular solution of the 
homogeneous equation 

i6M = [Ho, 6M] (52) 

to the above-obtained solutions. This can be achieved by turning the phases e^*"^* into 
gztiojt _ g-j(_Ej(k)-E^ (k ))t ^^ ^Yie above expressions ( |5DD and (^TJ). The solutions are explicitly 
given by means of the following formulae. For the matrix elements between two positive 
energy states t=+, r'=+ one has 

For the matrix elements between one positive and one negative energy state t=+, r' = — 
one has 

^g±J(a;T(Sa(k)+£;6(k')))i/2 _ (-54^ 

For the matrix elements between one negative and one positive energy state r = — , r' = + 
one has 



11 



Finally for the matrix elements between two negative energy states t = —^ r' = — one has 



.. ., -I.A..-.I., .', -)^ ,;g|,)!tl-)) ''-''''^-"''^'""'^' 



Here the indices I and m are summed over. At this stage, we are ready to give some com- 
ments. If temperature is equal to 0, then (a, fc, ib|(5A1|6, k\ ±) = because na{k) = 0. The 
dominant responses are {a,k,±\5Ai^^'^\b,k',^) because they have a two-meson pole. One 
can understand what physical process these terms represent. Let us remember the relation 
ttkln) = ^/nk\n — 1) and a[,|ra) = y/n^ + l\n + 1) for boson annihilation and creation op- 
erators and the eigenstates of boson number. The amplitude of two-meson creation minus 
two-meson absorption process is thus proportional to {nk + l){nk' + l)—nknk'. This is nothing 
but the factor appearing in (a, k, ±\SAi^^^\b, k', =i=). Thus a net two-meson emission is real- 
ized because they have the factor na{k) +ni,{k') + 1. Let us imagine the relativistic heavy ion 
or proton-proton collisions. The restored chiral symmetry spontaneously breaks down again 
after rapid cooling. The spontaneous symmetry break-down may then be realized at temper- 
ature T ^ with the explicit chiral symmetry breaking term. If so-called disoriented chiral 
condensed state is developed, this state should be finally relaxed to the normal condensed 
state due to meson emissions. In those situations, the response (a, k, ±\5Ai^^^\b, k', =p) will 
play an important role through two-meson emission from collective isospin rotating mean 
field (condensate) and will dominantly contribute to the damping time of the collective ro- 
tation in isospin space. We will return to this consideration when we estimate the damping 
time. 



VII. EFFECT OF INTERACTION TERMS 

In this section we take into account of the interaction term S7i which was neglected in 
the previous section. When this term is present the external field can induce a change ST 
in (|36| ) in the mean field given by the following formulae 

X N X 

S'^ab = ^^abYl "cc + -Oiab , (57) 

" c=l -J 

where the quantity a satisfies the following equation 

aab = {Xab + 5r„,)nlJ . (58) 
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Let us introduce the dressed response function Uik defined by 
for i = k wliile 



n..(^,q) = ^ f^\;;;^ . , (59) 



r(0) 






for i ^ k. With this notation we can write a in the following form: 



«22 - [X,, + Xu-Unj ■ u,, . ^_^n22-tnn' 

ail = (1^^ 2 in the above expression) , 

ai2 = Xi2ni2 , 

«21 = -^2in2i , 

aij = (for^,j>3). (61) 

The expressions for the physically important functions are thus obtained with a boundary 
condition 6A4 = at t = : 

(a, k, :f\SM^ '\b, k, ±) ^ , ' ,, , ^ ,, -- 

^ ' ' ^ lE,{k)E,{k±q) LO - {E,{k) + E,{k ± q)) 



x(xl?^+^rit))sm( -^^-(^)^^-(^^-^)) t) 

X 6{k' -kT q)e±^("+^"W+^*('^±'i))*/2 _ (g2) 

Here 5r^j,^(g) is related to the amplitudes X^^^ defining the external field via equations (|57D 
and (|61|). 

VIII. DAMPING TIME 

In this section we give the damping time of the collective rotating condensate in isospin 
space (or disoriented chiral condensation). We estimate the damping time by means of the 
ratio of the energy density of collective isospin-rotating condensate to the energy density of 
emitted mesons per unit time. To obtain the damping time, we need the energy density of 
meson fields originated from the "external source" term. 

The small change of the generalized density matrix is written by using a symplectic 
operator e"*-*^ as 

M = e'^^A^oe"'^ = Mo + SM + SM2 + ■■■ , (63) 

6M = z[x,Mo], 6M2 = -l[x,[x,Mo]]. 

Since we have known 6Ai in the previous section, we can get x as 
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{a,k,T\x\b,k',r') 



n{k)~n\k 



^r7r7z{a,k,T\5M\b,k\T'). 



(64) 



The variation of the energy comes from 6Ai2 because Trace(7io[X! -^o]) = due to 



AE = TTace(noSM2) 



akr fek'r 



V' ra{k)-n\k' 



(65) 



This energy shift is calculated by using the functions (a, A;, ±.\5M.^^'>\b^ k', =i=) because these 
terms mainly contribute to the energy shift satisfying the energy conservation as will be seen 
later. As a result we obtain 



AE, 



-^i<ekuu ^"(*) + «*') 



1 
ab - -- J ^ .^ J ^ '" na{k) + ribik') + 1 

x((a, k, +\6M^-^\b, k', -)(6, k', -\6M^+^\a, k, +) 
+ (a, k, -\6M^+^\b, k', +){b, k', +\6M^-^\a, k, -)), 



(66) 



where AE = J^ab^^ab- Here we can calculate the energy density by noting that the space 
volume V is expressed as V = J cPx = {2Tc)^6^{k = 0). 

The most important process is the two-pion emission corresponding to a = 6 = 2 among 
these responses. In this channel, we denote X22 = -^22 — -^22 and 6T22 = ST22 = 5T22- 
The energy density is expressed as 



AS^ 



AE^ 
V 



d^k 



n2{k) + n2{k - q) + 1] 



(27r)3 
y<iX22{q) + Sr22{q)) 



E2{k) E2{k - q) 



2 si n' [{u - E2{k) - E2{k - q))t/2] 
{to - E2{k) - E2{k - q)f 



(67) 



Using the relation (sin(u;t/2)/(cij/2))' k, 27ct6{u!) for large t, which ensures the energy con- 
servation, we find that the contribution of the two pion excitation to the energy density 
grows in proportion to time. 



AS^ 



i-K-K{q)t 



(68) 



where 



d^k 



/u rv, 
—--{n2{k)+n2{k-q) + l) 



+ 



E2{k) E2{k-q) 



x-(X22(g) + 6T22{q)f2'K6{uj - E2{k) - E2{k - q)) 



(69) 



gives the rate of the increase of the energy density associated with two pion excitation 
which carries net momentum q^ = {uj,q). Similar expressions 77rcr and 70-0- for the avr- and 
crcr-channels are also obtained. 



14 



The damping time r(g) of the collective rotating condensate in isospin space is simply 
estimated in terms of the energy density of rotating condensate and the energy density of 
emitted meson per unit time : 

r{q) = Sohiq) • (70) 

Here, the numerator is the energy density of the condensate originated from the collective 
rotation in isospin space : 

= l^lico' + q') , (71) 

where we have set the pion mass equal to zero in the lowest order approximation for r(g) 
according to the results obtained in Appendix A and we have neglected a correction from 
the quantum fluctuations which we had shown to be small (see section 8 of I). 
The denominator 

is the total rate of the change of the energy density due to two meson emissions. 



IX. DISCUSSION 

We estimate the associated damping time and the meson emission rate in this section. 
Let us consider a simple case with q = Q. The rate of the change of the energy density due 
to two-pion emission with three momentum q = is simply written as 



7^^(^,0) = — 
Stt 



X22(cU,0)h +^n22(cc^,o) 



Cj2 - 4M 



2 



X {2712 (I) + l) e{uj - 2M2) , (72) 

where we neglected nii(g) which corresponds to the aa excitation because two sigma mesons 
are too heavy in our consideration. For the crvr-channel at T = we obtain 



7^a(^,0) = — 
47r 

1 



|Xi2(a;,0)|M +^ni2(a;,0) 



2 



2kQe{uj - Ml - M2) , (73) 



ko = — V(^' - (^1 + ^2)2) (c^2 _ ^Mi _ M^ 



Let us give an order estimation of the damping time in the case of g = 0. We numerically 

estimate the typical damping time at a; = 2\/2M.„ where Mt,{^ 0) is the physical pion mass. 



The reason is that the factor uJ^ j JuS^ — 4M^ which appears in ( |70D with (^) has a minimum 
value if M^ is close to M^ due to the isospin-rotating effect. In this case, the energy density of 
collective isospin-rotating condensate is /^u;^/2 ^ (160MeV)^ where we used v^o = /tt = 93 
MeV (the pion decay constant). Note that the change in the potential energy density induced 
by the chiral symmetry breaking term is 2M^f^ where we used c = f-^M^. The validity of 
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our approximation scheme requires f'^u'^ /2 ^ 2M^f^, i.e., u 3> 2M^, which is reasonably 
satisfied for our choice of u. The damping time is roughly estimated as 

T = T^^^ 40 fm/c at cj = 2^/2M^ and q = , (74) 

where we used m^ k, M2 ^ 0, M^ = 140 MeV, m„ ^ Mi ^ M„ = 500 MeV and A = 86.7 
and we introduced the three momentum cutoff A = 1 GeV to calculate the response 
function 1122 • This is rather long compared with the collision time of heavy ion collision 
and/or the formation time of so-called disoriented chiral condensation about a few fm/c. 
When we take into account of the small finite g, the damping time becomes slightly longer. 

It is interesting to note that the calculation we have just made also provides information 
on the meson emission rate. Indeed, from equation (|74|) we find that the energy density 
decrease per unit time is of the order of (100 MeV)^ per fm/c. If we assume that the 
classical field configuration occupies a volume of the order of 1000 cubic fm this implies an 
emission rate of 100 mesons per fm/c. 

Further studies of the present model are worth being pursued. For instance it would 
be interesting to consider initial conditions for which there would be a coupling between 
collective rotations and small intrinsic oscillations. Comparing the damping times of rota- 
tions and vibrations separately would indeed teach us which of the two collective modes is 
selected by dissipation effects. Another interesting problem is to combine the effect of the 
expansion of the meson field, already considered in references and |l7|, with collective 
isospin rotations. 

ACKNOWLEDGMENTS 

Two of the authors (D. V. and Y. T.) would like to thank Prof. D. Boyanovsky and 
Prof. H. J. de Vega for valuable discussions and comments. One of us (D. V.) wishes to 
express his appreciation to the Japan Society for the Promotion of Science for support to 
visit the Yukawa Institute where this work was started. Y. T. is grateful to the members of 
the Laboratoire de Physique Theorique des Particules Elementaires (LPTPE) of Universite 
Pierre et Marie Curie for their warm hospitality. He also would like to express his sincere 
thanks to Nishina Memorial Foundation for financial support to visit the LPTPE. T. M. 
is grateful to the Laboratoire de Physique Theorique des Particules Elementaires (LPTPE) 
of Universite Pierre et Marie Curie for the support and the hospitality which he received 
during his visit. 

APPENDIX A: PIONIC EXCITATION IN THE CHIRAL LIMIT 

In this appendix, we calculate the pionic response to a weak external perturbation in 
the linearized approximation. It is shown here that we can obtain in this framework the 
pion mass which is consistent with the Goldstone theorem requiring that the pion mass is 
exactly equal to zero in the chiral limit. We introduce a source term Ja{x) acting on the 
second component of the A^-component scalar field with 0(N)-symmetry. In general we can 
extract the value of the square mass of the meson field from the functional 6{ip{x))/6J{y) 
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which is the inverse propagator in the hniit of the small source. Our task is thus to calculate 
the mean field value (y?) depending on J . 

The equation of motion for the mean field ip with a source term J is obtained as 



2 A \ ^ 

□ +m — — C/9 X (y9 (y9 = J, 



(Al) 



where m\^^ has the same form as @) by denoting (p a.s (p. If J = 0, then the above equation is 
reduced to (|l|) with (Q) in section 2. We denote the Feynman propagator S{x,x) appearing 
in rn?[x) as S^^\x,x) for the case J{x) = in this appendix. It is here noted that the 
propagator S"*^"^ is diagonal in the isospin quantum number. 
Let us take Jaix) to be of the form 



Ja{x) = e6a2 exp{iujt -iq- x) . 



(A2) 



The propagator S{x^x) is changed according to the change of m?{x)^ m?{x) -^ m?{x) + 
(5m^(a;), from the case J{x) = to J{x) ^ 0, namely, S{x,x) = S^^\x,x) + 5S{x,x). We 
write the change m^(a;) as 



6mlf^{x) = 5m^j, exp^iut — iq ■ x) . 
As a result, we obtain 5S{x, x) : 

6Sabix, x) = 6ml^, ■ 11;^°''' (g^) exp{iujt - iq ■ x) 



(A3) 



(A4) 



where 11)^^° (g^) is the polarization tensor which is expressed in terms of the Fourier transform 
of S^^\x,x) being S^'^\p) as 



For the source term 



we use the following Ansatz for 6m^ : 



(A5) 



Sm 



/O rj ■■■\ 
T] ■■ 
•■ 



(A6) 



An ansatz of the same form in the matrix structure is assumed for the quantity 6S{x,x). 
Then the linearized evolution equation for the small shift S(p{x) = p{x) — p^^^'^x) reads 

(-q^ + m^- ^(^(°) X if^A 6p + 6m' if^''^ - ^{S^ x <^(°) + (^(°) x 6ip)ip^'^ = J , (AT) 

where ip^^\x) is the value of the mean field in the case J{x) = 0. For the mass matrix we 
obtain 



Sm' = ^{5ipx (^(°) + ^(°) X Sip) + ^SSix,x) 



(A8) 
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which is consistent with the ansatz ( IKEj ) because ip'^^^x) = 5ai^o and the second component 
of 6ip{x) is only non-zero according to the source ( [A 21 ). It should be noted that for the source 
there is no change in Trace SS{x, x) and no change in (f'^{x) at least in the lowest order. 



By using the result 



and equation 



we obtain 



6m ip 



^,niO) 



A 2 



AttCpoI) 



6ip 



(A9) 



where we used ip^^^ ■ 5ip = 0. From the linearized evolution equation ( |^) and the above 
relation ( |A9| ) it is possible to express the change 6ip in terms of the condensate in the presence 
of the external source J. It is given by 



-q' + m' + ^^l ■ ^nr 



1 Att(poI) 

i - 3iil2 



6{p = J 



(AlO) 



The left-hand side is related to the dressed propagator. This provides a prescription to 
calculate the pion mass. Indeed from the above equation (|A10|) we obtain the J-dependent 
part of the mean field v?(x), that is, the shift Sip{x). This gives the inverse propagator for 
the meson field which is precisely the functional derivative 6ip{x)/6J{y). The square of the 
pion mass M^ is therefore just the value of the bracket calculated at g = 0, i.e., 



Mi 



/i 



^^l ■ ^nr^ 



Att(poI) ' 



• ■'^J-12 



(All) 



[SfUeansnSf)(g2 
following matrix structure [jl8[ : 



where 11^2"''' means 11^2" l*?^ = 0)- Here, we have used the fact that the mass matrix has the 



m 







V : 





















/ 



where the masses M and fi satisfy the self-consistency conditions : 

M^ = ml + ^^l + ^G{M-) + {N- 1)^G(/.^) , 

l^' = ml + -^^l + \g{M^) + (AT + 1)^G(^^) 

as well as the condensate minimization equation 

A 



M' 



-,vl 



(A12) 



(A13) 



(A14) 



The function G{m'^) is defined by 

G{m^) = 



(27r)V 



771^ + ie 



(A15) 



One can further use at this point the expression of the polarization tensor ( |A^ ) in terms of 
G{m^), that is, U'f^°^\q = 0) = [G{M^) - G{fi^)]/{M'^ - /x^). Also, fi"^ can be expressed in 
terms of G{im?). Namely, by subtracting the second equation from the first one in (|A13|) and 
by using the minimization equation ( |A14| ), we obtain jj? = — A[G(M^) — G(/i^)]/3. Finally, 



from ( |A11|) , we can derive the following relation by using the above-derived expressions of 
the polarization tensor and /z^ : 

(l-^nif^(g = 0))A42 = 0. (A16) 

This gives M^ = which is consistent with the Goldstone theorem. 
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